Constraints on the fundamental string coupling from B-mode experiments 
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We study signatures of cosmic superstring networks containing strings of multiple tensions and 
Y-junctions, on the cosmic microwave background (CMB) temperature and polarisation spectra. 
Focusing on the crucial role of the string coupling constant g s , we show that the number density 
and energy density of the scaling network are dominated by different types of string in the g 3 ~ 1 
and g s <C 1 limits. This can lead to an observable shift in the position of the B-mode peak — a 
distinct signal leading to a direct constraint on g s . We forecast the joint bounds on g s and the 
fundamental string tension pf from upcoming and future CMB polarisation experiments, as well as 
the signal to noise in detecting the difference between B-mode signals in the limiting cases of large 
and small g s . We show that such a detectable shift is within reach of planned experiments. 
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The extraordinary wealth of cosmological data of in- 
creasing precision is allowing us to seriously investigate 
models of the very early universe. Here we focus on as- 
pects of one promising model, namely the potential cos- 
mic microwave background (CMB) signatures of cosmic 
superstrings - small tension line-like objects stretched 
over cosmological scales - whose existence is characteris- 
tic of a class of inflationary models based on string the- 
ory. The observational properties of these strings depend 
sensitively on the value of the string coupling g s , a fun- 
damental parameter of string theory, and here we show 
for the first time that it may be possible to constrain this 
parameter through future CMB polarisation data. 

Cosmic superstrings differ from their topological coun- 
terparts, i.e. field theory cosmic strings in a number 
of ways. First, they have a spectrum of tensions fM de- 
termined by g s and the (pi, ft) charges they carry @, 
Hi = (hf 1 9s) \/ Pi 9s + 9ii where \ip is the tension of the 
lightest (1,0) fundamental string. Second, their inter- 
commutation probabilities can be significantly less than 
unity 4 6], so when two cosmic superstrings collide they 
do not necessarily reconnect. Third, when strings of dif- 
ferent charges collide, Y-junctions - i.e. vertices joining 
3 string segments of different charge - can form Q. In 
order to generalise previous work on the B-mode CMB 
polarisation signal from standard cosmic strings to 
cosmic superstrings, a first step is therefore to solve for 
the dynamics of a multi-tension string network with junc- 
tions. 

To do so we use the network evolution model developed 
in in which the fundamental variables are pi and Vj, 

respectively the energy density and root-mean-squared 
(rms) velocity of strings of type i, with correlation length 
defined by pi = PiL~ 2 . They obey the coupled equations 
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where • = d/dt and a(t) is the Friedmann-Lcmaitre- 
Robertson- Walker (FRW) universe scale factor. 

The presence of Y-junctions is reflected in the last two 
terms in (|2[); the penultimate one, for example, models 
the energy loss from the strings of type i due to their colli- 
sion with strings of type a (with average relative velocity 
Via), resulting in the formation of type k links of average 
length i\ a (t). This occurs with p s -dependent probabil- 
ity dJ- a = d^ which depends both on the microphysical 
quantum intercommutation probability Via calculated in 
and on a purely classical kinematic constraint Sf a 
described in @,[l(|. Self-interactions of strings of type i, 
leading to the formation of loops, are quantified by the 
coefficients a in © . These are determined by the micro- 
physical self-intercommutation probabilities Vu and are 
also <7 s -dependent. Finally, all the Vij also incorporate 
a model-dependent volume factor described by another 
parameter w pd| . which in this Letter we set to unity. 

Following the results [l2j of simulations of single-string 
networks with varying intercommuting probabilities (and 
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no junctions d\j — 0), we have taken a = c 
Matching to single-string network simulations with Va = 
1, then fixes c = 0.23 (0.18) in the radiation (matter) era 
[13| . There are at present no network simulations with 
which to calibrate interactions between strings of differ- 
ent types (i ^ j) , but one expects a similar dependence 



on Vij so we take d^ 
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Despite the infinite hierarchy of Eqs. (HHU, in practice 
only the first few lightest strings dominate [1(1 EH], 
so we truncate the system at i = 7 keeping all coprimc 
charge pairs up to (1,3) and (3, 1). The generic attractor 
solution of this system, with Vij and as computed in 
[Til ] . is a scaling one in which £j = Li/t and Vi asymp- 
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totically approach constant values. For all studied val- 
ues of the string coupling g s , the three lightest strings, 
namely (1,0) F strings, (0,1) D strings and (1,1) FD 
strings, dominate the string number and energy density 
of the network. This is due to the interactions between 
the strings, and is expected from microphysical studies 
of junction dynamics Q- Interestingly, however, while 
the hierarchy in the string number density JVj = re- 
mains the same (i.e. F strings are the most populous for 
all g s , followed by D and FD strings), the power spectrum 
density (llj of the network Mj = (l^i/^i) 2 becomes dom- 
inated by the heavy D strings at small string couplings, 
despite the fact that they are rare. In Fig. Q] we plot 
the rms velocity, the correlation length and the power 
spectrum density at the time of last scattering (LS) as 
a function of the string coupling g s for the three light- 
est network components. The power spectrum density 
is dominated by the light populous F strings for g s < 1 
and by the heavy rare D strings for j s < 1. This gen- 
eral trend appears to be quite robust [ll|, and is one of 
the main reasons behind the observable change in the B- 
mode spectrum. The other key ingredient is the fact that 
the correlation length of the D strings increases with de- 
creasing g Sl as shown in Fig.[T] This increase is greater for 
the weaker dependence of c& coefficients on Vij . Further 
investigation is needed to confirm the functional depen- 
dence of the dfj coefficients on Vij, and to understand 
better the non-perturbative interactions that determine 
Vij for heavy strings - the two main sources of uncer- 
tainty in our modelling. 

With the scaling solutions for interacting F-D networks 
at hand, we can predict their temperature (TT) and B- 
mode polarisation (BB) spectra. This is done using a 
generalised version of CMBACT [l5[ for multi-tension 
string networks modelled by equations ((HU), and im- 
posing that the total TT contribution from strings can- 
not exceed 10% of the observed CMB anisotropy [16h18I |. 
The amplitude of the CMB angular spectra Ci for these 
multi-tension networks can be written as the sum of the 
power spectrum densities Mi of each string type: 
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so this ultimately constrains the fundamental string ten- 
sion jjip. Here, we adjust fip to be such that 



f s = C. 
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where Curings an d ^totai are defined by summing the 
corresponding C/s over multipoles 2 < I < 2000 0]. 

Because of the important vector mode contribution 
from strings, the string induced B-mode signal can be 
much stronger than the one from inflationary tensor 
modes [1, and, indeed, strings can be a prominent 
source of B-mode polarisation on sub degree scales even 



with a marginal TT contribution of 1% [ll| . The shapes 
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FIG. 1: Dependence of the rms velocity (top), correlation 
length (middle) and power spectrum density (bottom panel) 
at the time of last scattering (LS) on the value of the string 
coupling g 3 for the F strings (solid black lines), D strings (blue 
dashed lines) and FD strings (red dotted lines). The rms ve- 
locity and correlation length of the string type(s) dominating 
the power spectrum are shown with oversized dots. 



network, namely the correlation length and the rms ve- 
locity. In particular, the position of the main B-mode 
peak moves towards lower I for larger correlation length, 
while the dependence on the rms velocity is nonlinear. 




of the spectra depend on the large-scale properties of the 
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FIG. 2: The B-mode power spectra for g s = 0.04 (solid lines) 
and g s = 0.9 (dashed lines) normalised to f 3 = 0.1. The insert 
plot shows the position of the peak as a function of g 3 . 

Fig. [2] shows the B-mode spectra for our two extreme 
values of the string coupling, namely g s = 0.04 and 
g s = 0.9. Our main result is a significant shift in the 
peak towards low multipolc number as the string cou- 
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pling is reduced. This has a physical interpretation in 
terms of the network behaviour. As mentioned above, 
the BB peak location depends on the correlation length 
and rms velocity of the string type dominating the power 
spectrum density Mj. Our results show (Fig. [1| that the 
velocity of the dominant string type only varies between 
0.6 and 0.7, so we expect the peak position to be de- 
termined mostly by the dominant correlation length £j. 
As we saw, the power spectrum density is dominated by 
F strings for large g s , but at lower g s there is a tran- 
sition and the rare (larger heavy D strings become 
more important. As we move to even smaller g s , the D- 
string correlation length increases (Fig. [TJ middle panel) 
and the peak moves further down in I. Our result im- 
plies that an observation of the B-mode peak could help 
to break the degeneracy between [ip and g s in the con- 
straint imposed by the normalisation condition (|4]) . This 
degeneracy may also be broken [III by combining CMB 
with pulsar timing constraints (20j. 
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FIG. 4: The S/N of an overall detection (lower panel) as 
a function of f a for a reference model with g s = 0.1. La- 
bels are the same as in Fig. [3l with the sets of three curves 
(showing the optimistic, realistic and pessimistic lensing sce- 
narios) from left to right. The upper panel shows the S/N of 
the difference signal between two models with g s — 0.04 and 
g s = 0.9, with these both having the same overall detection 
S/N of the reference model at fixed f B . 



FIG. 3: The dot-dashed (black) curves show an f s of (from 
top to bottom) 0.1, 0.01, 0.001 and 0.0001 in the (fi F , g s ) 
parameter plane. Other curves show the S/N = 5 contours 
of an overall detection of the B-mode signal for various ex- 
periments. The three curves for each experiment show the 
pessimistic (long-dashed lines), realistic (short-dashed lines) 
and optimistic (solid lines) scenarios for the removal of lensing 
B-modes, as discussed in the text. We also show the pulsar 
bounds based on data in [2(|] (dot-dot-dot-dashed lines) in the 
small (upper curve) and large (lower curve) string loop limits 
considered in [Til. [l8|. 

What are the chances of detection of the peak? We 
can estimate the ability of upcoming CMB polarisation 
experiments (using the relevant instrument sensitivities 
for each experiment) — Planck [H [H[, SPIDER [Hj], 
EBEX 1241 . Polarbear J25[ (using the parameters given 
in Ref. |26|]), QUIET |27| (using the phase 2 specifica- 
tion), as well as the proposed COrE [28[ mission — to 
detect the position of the peak in the string sourced BB 
spectrum. We also refer the reader to a recent work in 
j29j which studied the ability of the proposed CMBPol 
experiment (30j to distinguish between B-modes sourced 



by textures, Abelian Higgs cosmic strings and gravita- 
tional waves sourced by inflation. We follow the proce- 
dure for obtaining the combined noise power spectrum 
for each experimental channel from Ref. (30|. In Fig. [3] 
we show contours of constant signal-to-noise (S/N = 5) in 
the overall detection of the B-mode signal by each of the 
above experiments, depending on the values of the pa- 
rameters fi f and g s . Note that constraints on f s will be 
significantly tightened by measurements of TT, TE and 
EE from Planck and other experiments, but they cannot 
break the degeneracy between fip and g s . For reference, 
we plot contours of f s = 0.1, 0.01, 0.001 and 0.0001. 
Unfortunately the string B-mode signal can be obscured 
by those produced by gravitational lensing of the pri- 
mary CMB. The sensitivity to strings will therefore be 
affected by the extent to which lensing can be removed. 
We account for this in three scenarios - in the pessimistic 
case, we add the lensing component to the instrument 
noise power spectrum to obtain an effective noise, i.e. 
Np S = Ni + C l e cns . The lensing spectrum is computed 
using the CMB code CAMB for the best-fit WMAP7 
ACDM model 311. In the second, we assume the added 



lensing signal can be reduced by a factor of 7 (although 
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it may be more) in the white noise regime, as suggested 
by a quadratic estimator of the projected lensing poten- 
tial [321. Finally, we also show the most optimistic case 
where all the lensing is subtracted. We note that con- 
tours of constant S /N closely match contours of constant 
/ s . In the lower panel of Fig. 0] we therefore present the 
S/N as a function of / s for a fixed reference model with 
g s =0.1. This can easily be interpreted in terms of gen- 
eral up and g s by referring to Fig. [3] Note that even 
without lensing subtraction, B-modes from strings can 
be detected by QUIET at high S/N even if strings con- 
tribute to TT just above the 0.1% level. On the other 
hand, the string contribution would have to be close to 
the upper bound allowed by WMAP (f s ~ 0.05 - 0.1) 
to be detectable by Planck and SPIDER. The higher 
sensitivity of future experiments allows the possibility of 
distinguishing between models with different string cou- 
plings. We illustrate this by computing the difference 
signal between our two extreme values, g s = 0.04 (model 
A) and g s = 0.9 (model B). To determine the S/N, the 
tension fip is adjusted for each of the two models, un- 
til they have the same overall detection S/N dotoct of B- 
modes as the reference model (the latter obtained from 
the lower panel of Fig. [4j. The difference S/N is then 

computed via (S/N diff ) 2 = £ £ (Cf - Cff / (ACf B ) 2 , 
where the noise, including cosmic variance, is ACf B = 
y/2/(2£+l) (Cf + Cf + Nf) . We show this in the up- 
per panel of Fig. 01 Since the difference signal is al- 
ways more difficult to detect than an overall detection, 
S/N diff < S/N dctcct , with the degradation factor depend- 
ing on the experiment and lensing removal. We note that 



B-mode spectra are practically independent of other cos- 
mological parameters, especially within the range cur- 
rently allowed by WMAP, because B-modes are sourced 
primarily by vector modes generated during the brief 
window of the last scattering surface. 

We see that Planck and SPIDER arc unable to dis- 
tinguish between the two models at high S/N, even if 
the string fraction is close to the WMAP limit. There is 
a bigger degradation in S/N^s for SPIDER, due to the 
lower angular resolution (the smallest SPIDER beam has 
an FWHM of ~ 20', compared to ~ 5' for Planck), hence 
making it more difficult to observe a shift in the B-mode 
peak. For EBEX the two models could be distinguished 
at low significance if the string fraction is close to the 
WMAP limit. For the Polarbcar and QUIET experi- 
ments, it is possible to distinguish between the two mod- 
els with a > 1% contribution to TT, depending on the 
delensing efficiency. COrE improves this discovery space 
by roughly an order of magnitude, and hence will be 
a powerful observatory for testing fundamental physics. 
We conclude with the optimistic thought that within the 
next generation of B-mode experiments we may well have 
the means to pin down the value of the string coupling 
constant, a direct observational test of string theory. 
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